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Abstract 

We present a systematic expansion of Kramers equation in the high friction Hmit. 
The latter is expanded within an operator continued fraction scheme. The relevant 
operators include both temporal and spatial derivatives and a covariant derivative 
or gauge like operator associated to the potential energy. Trivially, the first order 
term yields the Smoluchowsky equation. The second order term is readily obtained, 
known as the corrected Smoluchowsky equation. Further terms are computed in 
compact and straightforward fashion. As an application, the nonequilibrium ther- 
modynamics and hydrodynamical schemes for the one dimensional Brownian motion 
is presented. 
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I. Introduction 

Kramers equation |I[] (although originally proposed by Klein [0]) and the associated 
Smoluchowsky equation [^ have been widely applied to diverse problems, namely: Brow- 
nian motion in potential wells, chemical reactions rate theory, nuclear dynamics, general 
stochastic processes, just to mention a few. A relevant overview (by no means complete) 
can be found in the review articles @]-01; books ^-0] and references therein. The pas- 
sage from Kramers to Smoluchowsky equation was informal and qualitative in the early 
days (see for example |1| and ||^). Only fairly recently a well founded relationship was 
established |T2[-P^. A multiple time scale approach has been put forward more recently 



2|,|4 



In section II we present a differential recursive expansion for Kramers equation, 
with respect to the inverse of the friction coefficient, within a continued fraction formal- 
ism p5[- p^. In a compact and straightforward fashion we retrieve the original and the 



corrected Smoluchowsky equation [^,|10|, in first and third order respectively. Higher 
order corrections are eased out with reasonable algebra, amenable to a diagrammatic 
framework. In section III some applications are presented, namely the hydrodynami- 
cal and non-equilibrium thermodynamical picture associated to the Brownian motion. 
Finally, in section IV concluding remarks include further applications and future work. 

II. Continued fraction expansion 

Kramers [|1| considers the probability distribution P(x, v,t) in phase space for the 
Brownian particle (mass m) in an external potential V^(x,t). The friction coefficient is 7, 
D the diffusion coefficient and T the bath temperature. Kramers equation (in the one 
dimensional case to simplify our analysis) is given by 

dt dx m dx dv J dv \ dv ) 



If / denotes a characteristic length scale (associated to V for instance) and with 
Einstein relation D = (ksT p^ {( is the mobility coefficient, (^^ = rri'-f), we introduce 
dimensionless variables p4| 

V{x,t) 



t ^ jkeT 



m 



X ^ m A, \ 
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and a dimensionless friction coefficient inverse e given by 



e = - ^ 
7/ V rn 

Thus, Kramers equation is cast as 
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Following |T^ we use the ansatz 



P(g,e,r) = Cexp('-M fl ^n(g,r)«l>„(0 



(3) 



where the $„(^)'s are the orthonormal Hermite functions [^ with \E'„(g, r) = $n(0 — 
for n < and C is an overall constant. Substitution yields a set of differential recursive 
(difference) equations for the \E'„s. Introducing the compact notation 



f^ ^ 
^=^' 



V = ^ Z^ = V + (VA) 

oq 



we obtain the recursion relation 29 



{ed + n) ^n{q, r) = -£v^D^„_i(g, r) - eVn + lV^n+i(g, r) (4) 

This recursion relation can be readily generalized to higher dimensions via the tensor 
Hermite polynomials |30]. The first relation obtained from equation (^) with ra = is a 
continuity equation, defined here as the generalized Smoluchowsky equation. 



9^o(g,r) + V^i(g,r) = 



(5) 



For the \l/n(a;, r) we may define an associated function \I/n(a;, s) via tlie transforma- 



tion 



*n(?,i") = / rfs^„(g,s)exp(-sr) 



(6) 



Both the integration path F and particulars of the s variable are to be defined 
adequately (Laplace or Fourier like), assuring among other things that the expansion 



n \ n / -i— ^ \ n / 
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is a well defined expression. Here we do not require the transformed functions, as in 



previous works |T3|-|2^. Furthermore we define 



ar. 



^/^DGn, Pn = Vn + IVG, 



yielding a compact form for the recursion relation 



^n(g, r) = a„^„„i(g, r) + ^„^„+i(g, r) 



With standard recursive relations and continued fractions techniques (p5|-|2^) the 
above differential recursion relation can be solved as 



^„(g,r) = £„(5,V,Z^)^„_i(g,r) n>0 
d^o{q,T) = -V^i(g,r) = -V£i^o(g,r) 



The last equation is an infinite continued fraction differential operator 



Cn{d, V, £>) = 1 - a„ ( 1 - an+l (l - a„+2 (1 ) ^ Pn+3) /?n+2 ) Pn+l A 
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This equation can be expanded into a sum of separable products of the kind 



Cn{d,V,D) = Y,BnAV,D)Rr,,i{d), Rn,m = J^ C^y'^^G^, ■ ■ ■ G^^ (9) 

1=1 ai,...ai 

On the other hand, from equations (||), (|^ and (|]) we obtain the identity 

G«^o = --(i--v/:T'vi/o (10) 

a \ a J 
Let us write the ansatz (an e expansion in purely spatial operators) 

oo 
Ln{V,D)=Y,ZnA^^^V (11) 



l=\ 



With this ansatz for ri = 1 we compute ([T0|) , then (^, obtaining an iterative relation 
for the Zi ;. Once the latter are computed we are able to compute all the £„s in terms of 
purely spatial operators. We present the expansion results up to c' for the first five \E'„s 

where 

Tg = V^A, Tg = 2 (VM)^ + ivM + (y^A)VA + -(V^A) V 

The generalized Smoluchowsky equation, from equation (|]) and up to e^ terms is 
therefore 

d^o =eV{l + e^Ts + e^T,) D^o (12) 



This expression coincides with previous results, see for example |]TB|. We believe 
to have presented a simpler and compact version for the recursion relations, amenable 
to further diagrammatic manipulations. Little attention has hitherto been paid to the 
evolution of moments higher than \i/o. Within an effective medium scheme [^ we assume 
\E'o to be known via the solution of (|T2D and thus compute the higher moments \E'„ (n > 0) 
in terms of \E'o- Iii this fashion we explore both formal and practical consequences (the 
usual scheme is to solve the recursion relations in a hierarchical form, \E'n in terms of 



III. Applications 

We show that the recursive relations (^) correspond to hydrodynamic balance equa- 
tions derived from the Kramers kinetic equation (for equations other than Kramers see 



3^-[Q). We may picture the Brownian motion problem as a one dimensional fluid. 
Restoring to the original (dimensional) variables, let us write some macroscopic densi- 
ties (moments) of interest |^, namely: the mass p (particle n), the momentum flow Jm 
(velocity density J and mean velocity u) , the hydrodynamic pressure 11, the kinetic energy 
E and the heat flow Je, given respectively by (hereafter we denote P(a;, f , t) = P{q, C,, r), 
given by (|)) 



p{x,t) = m dvP{x,v,t) = mn{x,t) 

J^^dvvP{x,v,t) 



p{x,t)u{x,t), u{x,t) 



roc 
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JE{x,t) = -m I dvv^P{x,v,t) 



The overall constant C in equation (^) can be adjusted such that we have the 
correspondence p{x,t) = '$o{x,t) (P, yielding C^^ = 2^/^T'mkBT. From the properties of 
the Hermite functions ||29ll we have 



p{x,t) = \E'o(x, t) = mn(x, t) 

Jm{x,t) = VT^i{x,t) = 'mn{x^t)u{x,t) 

E{x,t) = -U{x,t) = -vl(^^o{x,t) + V2^2{x,t) 

JE{x,t) = i=4Q^i(x,t) + v^^3(a:,t)j 

From the recursion relation (^, we obtain 
dp dJrn dpu 

dJm 



dt 
dE 

'dt 



dx dx 


K{x) = ^— 

m ox 


-^ + KJ^-2^{E- 


- Eq) , Eq = -v^p 



These are balance equations for the mass, momentum and kinetic energy density respec- 
tively. The RHS of these equations are the diffusive (gradient), convective (proportional 
to the external force) and dissipative terms, respectively. The first equation is the conti- 
nuity (generalized Smoluchowsky) equation, the only one having a conserving character 
(no dissipation). There is a convective (drift) term hidden in \E'i. The dissipative terms 
appear in a relaxation time form. The inverse relaxation times are 0,7,27, ...727 (as we 
compute higher order terms). Here we connect to the recent multiple time scale analysis 



mentioned in the introduction |^,||2^. Our formalism can be readily generalized for 
multiple components undergoing chemical reactions, provided the reaction terms are cast 
in the linear approximation. This is the case for inhomogeneous semiconductors in the 
semiclassical approximation |3H]-|3f 



We write down the Boltzmann entropy density (hereafter A;^ = 1) 
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a(x,t) = — dvP(x,v,t)\n — P(x,v,t) = dva( 
J me J 



x,v,t) 

with In e = 1 . This entropy density is shghtly different to the interpolative entropy defined 
in [^]. The equihbrium state is defined as the stationary state 9\l/^^ = 0, n = 0, 1,2 ■ ■■, 
where all flows have ceased, i.e. \1/^^ = 0, n = 1, 2, ■ ■ ■. In this context the equilibrium 
state is properly defined by D'^'^ = 0, and we readily obtain (see |^ for a one dimensional 
ideal gas) 



(TcqKq(a;),T) = neq(x) f 2 ~ ^"^ ^w^j j ^Q*^^) = S~h?~ 

Some nonequilibrium stationary states (ss) may exist and will be reported elsewhere. 

With the above defined entropy, our hydrodynamical picture allow us define the local 
nonequilibrium densities for the energy E., free energies F and G; the kinetic temperature 
9 and the internal chemical potential /lint as 



E{x,t) = -e{x,t)n{x,t) 

F{x,t) = E{x,t)-e{x,t)a{x,t) 

G{x,t) = F{x,t) + Il{x,t) = fii^t{x,t)n{x,t) 



These definitions do not assume the local equilibrium assumption |^, |3^ and will 
prove to be consistent. The total chemical potential is fi = fii^t + V{x) |^. From our 
hydrodynamical results, and hereafter to second order in e, we obtain for the kinetic 
temperature and the particle flow density 



^ I nD'^n(x,t) 

e{x,t) = T[l + e^ ^ ' ' 



n{x,t) 
J{x,t) = —evTDn{x,t) 



We estimate the temperature fluctuations \9 — T\/T for constant external fields, for 
two cases at room temperature. First, for an Aluminum Brownian particle (radius ~ 10^^ 
mm) in water, the constant force being gravity. The temperature fluctuation is of the 
order of 10~^. Second, for a typical semiconductor under an external electric field [^ ~ 10 
Vcm~^, temperature fluctuations are of the order 10~^. In the latter case, temperature 
fluctuations reach a few percents of the bath temperature only for extremely large fields 
of the order of the breakdown field ~ 10^ Vcm^^ where the semiclassical picture may fail 
for a typical semiconductor. 

The entropy is given by 



a{x,t) = creq(n(x,t),6'(x,t)) - 
= aieq{x,t) +ag,{x,t) 



2n{x,t) 



d_ 
dx 



+ A{x,t) 1 n{x,t) 



where we have separated the local equilibrium contribution from the gradient correc- 
tions term, the latter with a Ginzburg-Landau reminiscence ^^ , here with the covariant 
derivative (the external force is the affine connection). Indeed, the local equilibrium ap- 
proximation is valid up to linear corrections in e. The total chemical potential is given 
by 



fj,{x,t) 

ip{x,t) 



^{x,t) + 
e{x,t)\n 



e^ e{x,t) 

Tn2(x,t) 

n*{x,t) 



d 

— + A{x,t) ]n{x,t) 



n*{x, t) = n{x, t) exp 



V{x) 

e{x,t) 



As expected for the total chemical potential p8| 

dfi{x,t) 



J{x,t) = —(n{x,t)- 



dx 



Finally we multiply (|T]) by the integrating factor 1 + In (;j^P(x, f , t)) and integrate 
over the velocity. We obtain the entropy balance equation (integrating by parts, we 
assume v"'P{x,v,t) ^ 0, v —^ ±oo). 



daix,t) dJ„{x,t) _ 

dt ^ dx -^-^"^'^^ 

The computed entropy flow density J„{x,t) (associated to entropy exchange with 
the bath) and entropy production Pa{x, t) (associated to internal entropy production) are 
given respectively by 

Jfj{x,t) = dvva{x,v,t) = a{x,t)J{x,t) 



Pa{x,t) = -f dvP{x,v,t) [l-jl 



din P{x,v,t) 



2' 



dv 
= ^ {E{x, t) - Eo{x, t)) + 27 {a{x, t) - aieq(x, t)) 

The entropy production density is related to both an energy and entropy gradient 
'excess' in a relaxation time fashion. 

From Smoluchowsky's equation (|1^) we also have the identity 

^, , dnixA) ^, .dhmixA) 



IV. Concluding remarks 

The recursive method presented is compact and amenable to diagrammatic tech- 
niques (work in progress, as an e expansion). We did show that the moment (\E') recursion 
relations are nothing but the hydrodynamical balance equations. Therefore we can read- 
ily obtain the diffusive, convective and dissipative contributions for any relevant quantity 
such as particle number, momentum and energy density. 

Furthermore, we presented a nonequilibrium thermodynamical picture of Brownian 
motion, hinging on the knowledge of the fundamental moment, namely the density p, 
calculated via the generalized Smoluchowsky equation. 

9 



With the introduction of a generahzed entropy we derived the nonequihbrium ther- 
modynamic density functions, the kinetic temperature and the chemical potential. We 
obtain gradient corrections to the local equilibrium approximation. We presented some 
relations among the computed quantities, namely: particle flow to chemical potential and 
entropy gradient corrections to entropy production. 

Work in progress will probe and renormalize the thermodynamics presented here, 
as higher order terms in the e expansions are considered. Other interesting applications 



are related to nonisothermal oscillatory chemical reactions p9| , ^ 
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